Preparation and manipulation of a fault-tolerant superconducting qubit 
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We describe a qubit encoded in continuous quantum variables of an rf superconducting quantum 
interference device. Since the number of accessible states in the system is infinite, we may protect 
its two-dimensional subspace from small errors introduced by the interaction with the environment 
and during manipulations. We show how to prepare the fault-tolerant state and manipulate the 
system. The discussed operations suffice to perform quantum computation on the encoded state, 
syndrome extraction, and quantum error correction. We also comment on the physical sources of 
errors and possible imperfections while manipulating the system. 

PACS numbers: 85.25.Cp,03.67.Lx,03.67.Pp 



I. INTRODUCTION 

Maintaining quantum coherence is crucial for 
quantum-information processing. Coherence of any 
quantum system is gradually suppressed due to un- 
wanted interactions with the environment. Among 
proposed realization of qubits, solid-state devices appear 
particularly promising due to the scalability and ease 
of integration in electronic circuits, but their operation 
requires keeping them coherent, a potentially strong 
problem due to the host of microscopic modes. One 
of the achievements of quantum-information theory 
is the development of quantum-error-correction tech- 
niques, which allow one to suppress the effect of the 
environment on the software level, i.e. by running the 
appropriate quantum code to eliminate errors^ The 
standard approach protects against large errors that 
occur rarelyiSi 3 .^ Alternative methods were proposed 
by Gottesman and Kitaev — , where the shift-resistant 
codes protect against errors that occur continuously 
but are weak. Specifically, they analyzed the situation 
when a qubit is embedded in an infinite-dimensional 
Hilbert space of a physical device with a continuous 
degree of freedom. Continuous-variable quantum 
codes have been further developed for both qubits and 
qudits (d-dimensional analog of qubits) The 
continuous-variable codes became also a framework for 
discussion of quantum key distributionii or quantum 
teleportation of continuous quantum variables.—. 

In the original work^ the authors focused on encoding 
a qubit in an oscillator and developed codes which pro- 
tect against small shifts in the canonical variables of the 
oscillator. They also discussed implementation of this 
approach in optical systems. A universal set of quan- 
tum logic gates and error-correction steps may be realized 
using linear optical operations, squeezing, homodyne de- 
tection, photon counting, and nonlinear couplings. While 
some of these steps are easily realized with optical means, 
the encoding and certain gates from the universal set, 
which require photon counting and nonlinear couplings, 



are difficult to implement (alternative method of encoded 
state preparation has been discussed in Ref. [HI). 

Here, we suggest the implementation of the shift- 
resistant codes in superconducting devices. We show that 
their physical properties simplify the implementation of 
the difficult steps mentioned above. In these systems one 
may use the charg o 14 ' 15 ' 16 or the conjugate phase degree 



of freedo m 17 ' 18 to store and process quantum informa- 
tion. Typically, one adjusts the parameters such that 
the ground state is (almost) doubly degenerate, and at 
low energies the system reduces to a qubit. Here we con- 
sider a different approach, in which a continuous degree 
of freedom is used: the magnetic flux through the loop of 
an rf-SQUID (superconducting quantum interference de- 
vice) (a superconducting loop interrupted by a Josephson 
junction) and the conjugate charge. The Josephson cou- 
pling can be tuned if one replaces the junction with a dc- 
SQUIDfi^ For a turned-off Josephson coupling, the sys- 
tem behaves as a harmonic oscillator, and one can encode 
a qubit in its Hilbert space and manipulate the qubit 's 
state. The specific useful features of this design are re- 
lated to the periodic flux dependence of certain physical 
properties: this allows one to control the Josephson cou- 
pling and to monitor the magnetic flux modulo the flux 
quantum, thereby projecting out a comblike state needed 
for encoding of the qubit. 

The goal of our work is twofold. On one hand, we 
suggest an implementation of certain steps needed for 
shift-resistant codes, which are hard to implement with 
optical means. On the other hand, one may think of 
implications of our results for the long-term strategy of 
quantum computing in superconducting systems. While 
fabrication of larger circuits should be relatively straight- 
forward, it will be necessary to protect them against de- 
coherence. Although the approach we discuss here im- 
poses certain constraints on the system parameters and 
requires complicated operation procedures, the alterna- 
tive of standard quantum-error-correction codes requires 
building circuits with many auxiliary qubits and perform- 
ing complicated series of logic gates during error detec- 



2 



tion and recovery (cf. Ref. 20). Here, we demonstrate 
that the ideas of continuous-variable codes may be, in 
principle, implemented in superconducting circuits, un- 
derline advantages of these devices but do not optimize 
the design, and only superficially consider constraints on 
the circuit parameters. 

We begin with a short description of the shift-resistant 
codes in Sec. [TTJ In Sec. IIII1 we briefly describe the phys- 
ical system representing the qubit and discuss the pro- 
posed implementation of the difficult steps, i.e., of the 
encoding procedure and the quantum gates, including 
those needed for error correction. In Sec. IIV[ we discuss, 
for completeness, further necessary steps such as squeez- 
ing and translations in the phase space and two-qubit 
operations. We then comment on the error models rele- 
vant for practical devices and discuss constraints on the 
circuit parameters. 

II. SHIFT-RESISTANT CODES FOR LC 
OSCILLATOR 



For the purpose of error correction one measures the 
value of <i> mod 1. Such a measurement provides informa- 
tion on the possible error, a shift of the comblike states 
[Eq. JSJ], but does not distinguish between the combs 
peaked at even and odd flux values, |0) and |1). The 
observed shift is compensated by the inverse flux shift, 
which moves the comb structure to the closest integer. 
This procedure can correct sufficiently small flux-shift 
errors, A$ < 1/2. Similarly, one can correct small shifts, 
AQ < 1/2, in the charge variable Q. If the protection 
from rare large errors is also desired, these codes can be 
concatenated with the standard error-correction codes 4 

The states [Eq. ©] define a two-dimensional subspace 
protected from decoherence. In the following sections, we 
will show how one can prepare such states, store them 
(in the oscillatory regime with the Josephson coupling 
turned off), implement quantum logic gates and error- 
correction steps using the control over the Josephson 
coupling, the flux bias, and inductive coupling between 
qubits. 



While qubit is the simplest nontrivial quantum sys- 
tem, many physical systems offer the opportunity to uti- 
lize many levels and often a continuous spectrum to pro- 
cess quantum information. Gottesman and Kitaev — 
suggested to encode a logical qubit in a system with a 
continuous degree of freedom, an oscillator. They de- 
scribed error-correcting codes, which protect the state of 
the qubit against perturbations that cause weak diffusion 
of the position and momentum of the oscillator. Here, we 
briefly describe this proposal in the language of an LC 
circuit. 

Consider an LC oscillator, with the Hamiltonian 



$ 2 Q2 
H - 2L + 2C- 



(1) 



Here, the dynamical variables are the magnetic flux $ 
and the conjugate charge Q. For convenience, we use 
below the dimensionless flux and charge variables, 
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1/4 



Q: 



(2) 



and omit the primes. The rescaled variables satisfy the 
standard commutation relation 

[*,<?]=*■ 

In the proposal of Gottesman and Kitaev — the code- 
words (the states that encode the basis logic states of the 
qubit) are comblike superpositions, both in the flux and 
charge representations, 

oc oo 

io) = i<f = 2 S )= J2 \q = ™)> 

S — — QO Til — — OO 

OC oo 

|1) = J2 l* = 2«+l>= J2 (-l) m |Q = m).(3) 



III. ENCODING THE QUBIT AND LOGIC 
GATES 

A. Oscillator and the identity operation 




FIG. 1: The simplest flux qubits (Ref. \2w can be used to en- 
code a qubit in their continuous variables, (a) The rf-SQUID, 
a simple loop with a Josephson junction, and (b) rf-SQUID 
with tunable Josephson coupling. 

The Hamiltonian of an rf-SQUID with a tunable 
Josephson coupling [shown in Fig. [ljb)] reads 



H = TTUJo 



-£j($ x )cos 



91 

2 



2^(L/Cj)^§- 



(4) 



Here, the first two terms, which form the Hamiltonian of 
an oscillator, are the magnetic energy, controlled via the 
external flux bias <I> X and the inductance L of the loop, 
and the charging energy (depending on the parameters 
the capacitance Cj may reduce to that of the junction or 
involve the geometry of the whole loop), luq = (LCj)" 1 ^ 2 
is the frequency of the LC oscillator. The last term is 
the Josephson energy, with magnitude 



Ej($ x ) = 2E°j cos [ 7 r 3 / 2 (L/C J ) 1 / 4 $ x /$ c 
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controlled by the flux <I> X through the small loop of the 
dc-SQUID. 

We observe the qubit's states [Eq. ([3])] in the interac- 
tion picture of dynamics. In this approach, the state of 
the qubit is conserved if no manipulations are performed 
(although in the Schrodinger representation, the Hamil- 
tonian generates evolution), thus providing the identity 
operation. Specifically, consider the system in the oscilla- 
tory regime, when the Josephson coupling and the exter- 
nally applied flux $ x are tuned to zero. Then for an arbi- 
trary state, the Hamiltonian [Eq. |T])] , with the parabolic 
potential centered at the origin, generates rotation (of the 
density matrix in the Wigner form) of the state with the 
frequency uj in the phase space <&-Q around the origin. 
In other words, the amplitudes Qr in the rotating 
frame (interaction representation), which are related to 
<I> and Q as 



Qr 



COS LUot — shiuot 
sin cos wot 



(5) 



are integrals of motion, <&r = Qr = 0. Below we discuss 
quantum logic transformations in this reference frame. 
Note that performing operations in the rotating frame 
implies in practice that one has to keep track of the phase 
of the oscillator at frequency loq and perform all steps 
relative to this reference signal. 



B. Encoding 

Using the possibility to measure the flux mod <i>o, one 
can prepare a code state [Eq. ©], |0)or |1), by first 
preparing a state with a wide uniform distribution of flux 
and then measuring the flux value. As a result, one ob- 
tains a comblike state, with equidistant peaks and an 
offset of the structure, given by the result of the flux 
measurement. One performs a flux shift to compensate 
this offset and arrives at a codeword [Eq. Possible 
methods to implement these steps are discussed below. 

Preparation of a wide flux distribution may be real- 
ized by a number of methods. For instance, one may let 
the oscillator relax to the ground state in a narrow well, 
with the inductance decreased by (for instance, slowly) 
switching on a strong Josephson coupling. This creates a 
"squeezed" state, with a narrow flux distribution and a 
wide charge distribution. When the Josephson coupling 
is turned off, this distribution in the $-Q-plane starts 
rotating about the origin and after a quarter of the os- 
cillator period, ■k/{2u)q), a wide distribution of flux is 
reached, ~ J d$ |<f>). Alternatively, one can prepare such 
a state by squeezing the ground state of the oscillator 
(see below). 

Now, to single out a comblike structure out of this 
wide distribution, one is to perform a measurement of 
the flux value up to a multiple of a certain offset. This 
can be achieved by coupling the rf-SQUID inductively to 
a readout dc-SQUID and reading out the critical current 



of the latter. This critical current is a periodic function of 
the total magnetic flux in the measuring SQUID's loop, 



2I c0 cos 



where $ mx is the external flux bias of the measuring 
SQUID and the coupling A depends on the self- and mu- 
tual inductances of the loops, i.e., on the geometry of the 
SQUID's. For the particular geometry depicted in Fig. [21 
the parameter A can be close to 1 (our scheme is valid 
for different scenarios and values of A as well). 

Reading out the value of the critical current, one 
projects out a comblike state, 



E 



1$ = 



a. 



(6) 



with a period <!>p = $o/A (if only the absolute value of 
critical current is measured and the sign is not resolved) . 
a here is the initial displacement of the state from zero. 
(In fact, one would obtain a superposition of two comb- 
like states with different offset values a since two series 
of delta peaks correspond to each value of the critical 
current; one method to leave only one of these is to read 
out the critical current again after a small shift of the 
flux bias <E> mx .) 

Starting from the state projected by Eq. © from the 
wide flux distribution, one could compensate for the shift 
a by performing the inverse flux shift and also change the 
period of the comb structure by squeezing (see the next 
section) to arrive at one of the basis states, |5) or As 
we shall see later, for the purpose of performing one-qubit 
gates it is convenient to tune the peak separation to the 
superconducting flux quantum, $0 (possibly by squeez- 
ing). However, according to the definition of the encoded 
states [Eq. ((3}] and in the units defined in Eq. @ the 
peak spacing in the state |0) equals 2(7r 2 L/Cj) 1 / 4 . Com- 
parison to the flux quantum gives an extra constraint on 
the ratio L/C„ 
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L 

Cj 



1/4 



— = * /(2V5F). 



(7) 



With this condition satisfied, the Josephson term in the 
Hamiltonian of the system [Eq. Q] becomes periodic in 
$ with the period 2. 



C. Quantum gates 

If the Josephson coupling of the SQUID loop is 
nonzero, the Hamiltonian in the interaction picture has 
for $2, = the following form 

H = — i?j($ x ) cos {tt&r cos ujot + ttQr sin u; i] . (8) 

For the times t z — kir/uo and t x — ir(k + 1/2)/u)q (k 
being an integer), it reduces in the code subspace to 

H(t z ) = -Ej{^)a z , (9a) 
H(t x ) = -Ej($ x )a x . (9b) 
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FIG. 2: A possible method of the inductive meter-qubit cou- 
pling. In this geometry (when the qubit is placed inside the 
meter), the coupling is strong, and A ~ 1 . For simplicity, the 
qubit is shown without the smaller, dc-SQUID loop. 

Here, the Pauli matrices act on the encoded states |0) and 
|1) like on the spin states "up" and "down" along the z 
direction respectively. Using short pulses (the duration 
t should satisfy tloq <C 1) of the magnetic flux $ x at the 
times t x i z \, we can perform small phase shifts between 
the encoded states or induce transitions between them. 
Sufficiently large shifts, accumulated during many pulses, 
give rise to simplest one-qubit operations that generate 
the group of all possible unitary transformations in the 
encoded subspace^ 

Certainly no operation can be performed instanta- 
neously. So, apart from the instantaneous Hamiltonian 
corresponding exactly to the times t = t x i z \, we should 
analyze the time dependence of the Hamiltonian near 
tx(z) t° evaluate the error introduced by the finite rota- 
tion during the operation. Let us consider the Hamil- 
tonian [Eq. ^j] at times close to t z . Expansion to the 
quadratic terms in ujqt gives 

H w -Ej($ x ) {cos[7r$ fl ] -7rQnsin[7r$n]w T- 

-[(7rQ ii ) 2 cos(7r$j i ) - (ir$ R ) sin(vr$ jR )](cj r) 2 } . 

The error is generated by the time-dependent part. Since 
for the code space 

sin(7r$ji)|0) =sin(7r$ fl )|l) =0, 

the error is generated by at least second order terms in 
luqt, and for loqt <C 1 it is small. In addition, as shown in 
Sec. El errors generated by quadratic terms in Q can be 
corrected using the error-correcting routines. (The same 
property holds at times close to t x .) 

To complete the set of universal quantum gates, we 
need also a many-dimensional (in this case continuous) 
equivalent of CNOT, the SUM gate, which transforms 
the variables of two coupled qubits according to 

SUM: Qi->Qi-Q 2 , 

Here, the indices number the coupled qubits. This gate 
belongs to the so-called symplectic group and may be re- 
alized using (in quantum-optical setting) phase shifters, 
squeezing, and beam splitters^^ - elements that are ac- 
cessible also for LC oscillators, as we will show in the next 
section. The SUM gate reduces in the code subspace to 
CNOT. Its continuous nature is, however, crucial during 
the syndrome extraction. 



IV. MANIPULATING THE SYSTEM 

So far, we presented the computational steps that are 
difficult with quantum-optical elements and much more 
natural with the Josephson-junction systems. To make 
the discussion complete, let us now discuss the opera- 
tions that are necessary for the symplectic operations and 
can be performed on arbitrary wave function: squeezing, 
phase shifts, translations in the $-Q plane, and inductive 
coupling of two oscillators. In this section, all operations 
are described in the oscillatory regime Ej(& x ) = 0. 



A. Squeezing and phase shifts 

The scaling factors in Eq. ([2]) define dimcnsionlcss vari- 
ables for which the Hamiltonian is parametrized by only 
one real parameter, the energy scale luq. Since the scal- 
ing depends on the ratio L/C} and the frequency luq on 
the product of L and Cj, modification of either of the 
parameters (L or Cj) influences the system behavior in 
two ways - the oscillation period changes and the evo- 
lution (rotation) path in the phase space is squeezed in 
one direction. In particular, suppose that there is an ex- 
tra capacitance in parallel to the junction that can be 
switched instantaneously on and off, so that the total 
capacitance may equal Cj or XCj. Switching at some 
instant of time [for simplicity we assume that this mo- 
ment corresponds to t — in Eq. |(SJ)] from Cj to ACj 
rescales the variables, $ = — > A 1 / 4 ^ = A 1 / 4 ^, Q = 
Qr — * A _1 / 4 Q = A~ 1//4 Qi?., and modifies the frequency, 
luq — > u>i = A -1 / 2 ^. Using this property, we can perform 
squeezing of arbitrary wave function. However, according 
to Eq. ([7]), there is a constraint relating the system pa- 
rameters to the value of the flux quantum and we have to 
switch the capacitance back to the original value before 
performing further steps, like the single-qubit operations. 
The second switching (back from ACj to Cj) should be 
performed after a time delay, corresponding to an ex- 
change of the role of the variables. To be more specific 
suppose that we switch the capacitance from Cj to ACj 
at t = 0, switch it back at t = tt/2lui (quarter of the 
full oscillation), and observe the effect after the period 
of oscillation is completed (then the laboratory and ro- 
tating frames coincide and the effect is identical, which 
simplifies the analysis). The evolution in the Heisenberg 
picture of dynamics (in the laboratory frame) is described 
by the inverse of Eq. (0 , 

(S)=^W')(^). 

where U^ 1 is the inverse of the rotation matrix in Eq. ([5]). 
In this picture (and frame), the scaling of variables corre- 
sponding to the switching from Cj to ACj is described by 
a diagonal matrix S with the elements A 1 / 4 and A -1 / 4 . 
Thus, the entire procedure is described by the following 




& T squeeze = U 1 {uJq, 3ir/2uJ )S 1 U 1 ((J 1 , TT /2(J 1 )S 

~ ^ A" 1 / 2 ) ■ 
The interaction-picture result is 

This equation clearly describes squeezing of an arbitrary 
state. Note that after we change the capacitance of the 
qubit, the frequency is modified as well. This can be 
used in implementing another element - the shifter of the 
relative phase between two oscillators. If we do not want 
the state to be modified but only the evolution advanced 
(or delayed), we should switch the capacitance to XCj 
for the full period 2tt/ uj\. 



B. Translations in the phase space 

The evolution in the interaction picture (the rotating 
frame) is trivial, provided that the external flux <i> x is 
zero or at least constant. Once we switch it to a finite 
value, the physical center of rotation in the phase space 
shifts by $ x , and to simplify the description again we 
would have to shift the frame as well. However, square 
pulses of the flux <E> X give us the possibility to perform 
another useful and necessary in error-correcting routine 
transformation - translations in the phase space (in the 
<5-Q plane). 

To avoid these noninertial effects of the rotating frame, 
we describe this procedure in the laboratory frame (it 
lasts, however, for the full oscillation period, and like 
in the case of squeezing the effects are in both frames 
identical after the operation is over). Suppose that we 
begin with a state described by the wave function tpo ($) 
in the potential centered at the origin ($> x — 0). At the 
time t = 0, we turn on the external flux $> x = —(3/2. 
After the time t — tt/luq, which corresponds to a half 
of the oscillation, the state is transformed into ipi($>) — 
ipo(— $-1-/3). Then, we turn the external flux off and after 
the oscillation is completed, we obtain ^ 2 (^) = ipo(® — f3) 
- the initial state shifted by /3 in <I>. To shift the state in 
the Q direction we need to perform the same procedure 
on the Fourier-transformed state, i.e., to delay the initial 
moment of operation by tt/2ujo- 



C. Inductive coupling of oscillators 

Two qubits may be coupled inductively using an addi- 
tional LC circuit (see Fig. [3]) . The Hamiltonian of the 
LC oscillator and the qubits (for a moment we turn back 



FIG. 3: Inductive coupling of two qubits is realized with the 
Josephson energy tuned to zero. It is thus the coupling of 
two LC oscillators. In the rotating frame of reference, this 
coupling provides equivalent of quantum-optical beam splitter 
(see the text for explanation). 

to the natural units) is 

H = H 1 +H 2 + $ 2 /2L osc + Q 2 /2C osc -VQ, 

V = £m<*</Z, (10) 

i 

where Hi and H 2 are oscillatorlike Hamiltonians of the 
form of Eq. (JTJ), Q is the charge on the leads of the ca- 
pacitor C osc , and $ is the flux through the external LC 
circuit. If the frequency of the oscillator is much bigger 
than the frequency of the qubits ojlc = ]-/ \/L osc C osc ^> 
ujo = \j\fLC~i (in the limit C osc — > 0), we find that 
the coupling circuit remains in its ground state and the 
interaction (mediated by virtual excitations of the cou- 
pling circuit) is described by the term — C osc V 2 /2. Since 
$j = Qi/Cj, the interaction part of the Hamiltonian for 
two qubits has the form 

H mt = -J2 C °2^ j QiQj = -KizQxQi, (11) 

i=£j J 

so that after integrating out the oscillator's variables, the 
Hamiltonian of the interacting system equals 

H = H 1 +H 2 -K 12 Q 1 Q 2 . 

(Here, the parameters in i?i(2) are slightly different than 
in Eq. (|TU|) - from the term —C osc V 2 /2 we obtain the 
interaction part of the Hamiltonian (|11[) as well as terms 
quadratic in <3i(2), which formally slightly modifies the 
frequency of the qubits. In any case, this does not influ- 
ence the general shape of our results - the phase shifts 
between groups of qubits can be compensated using the 
previously discussed techniques.) We can rewrite this 
Hamiltonian in the rotating frame of reference, 

H = ^[QfliQfl2(l + cos2w f)+ 

$m$H2(l ~ cos2w i) + (Q m $ 2 + Qr2$i) sin 2uj Q t] . 

Since over time scales much longer than 1/cjq the effect 
of the oscillating terms averages out, we finally arrive at 

H = f- (QmQ R 2 + ■ 
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Now, we solve the Heisenberg equations of motion for 
the operators with the initial conditions $1(2) (0) = ^?( 2 ) > 
Qi(2)(0) = Q?(2) (f° r simplicity we omit here the R sub- 
scripts), and we arrive at 

/ » n #12* _n #12* 

^i(2)(t) = *; (2) cos + Q° 2(1) sin 



!l(2) 



(*) 



2 

Kiat 



- $ i(2) si n ., 



2(1)1 

V2(l) COS . 



If we make further the substitution (which is done in 
practice by advancing the second oscillator by 7r/2), 

$2 -> -Q2, <?2^$2, 

we obtain equations describing the action of a beam 
splitter— with time-dependent reflectivity and transmit- 
tance, 



#12* ■ #12* 

cos , yl = sin . 



R 



This element completes the set of all necessary opera- 
tions. 



V. PRACTICAL CONSIDERATIONS 

A. Error models and error recovery 

Errors that can be corrected using shift-resistant quan- 
tum codes are translations that leave the state out of the 
code subspace but still not too far from the initial state. 
The errors that can occur are due to the following inter- 
action with the environment: 



|a;) 9 b|0)_E — > (t 



(12) 



Certainly, almost no physical interaction is of this kind, 
but it can be expanded in terms of such translations, 



\x) qb \0) E -> / da'dp'C(a',0',t) (e iQa ' e 1 ^ \x)) \ct,p) E 

The syndrome extraction procedure gives then one pair 
of real parameters (a, (3) only, projecting the state onto 
Eq. ([12")) . Effectiveness of the error recovery depends in 
this case on the amplitude C . We may treat the function 
|C| 2 as the probability distribution of possible errors (or 
equivalently C as two-dimensional wave function). If the 
uncertainties Aa' and A/3' are small after the interval 
between two error-correcting routines, the states will be 
recovered with high fidelity. 

The interaction with the environment can be also writ- 
ten in the operator-sum representation (more convenient 
for our purposes) 

p{t) = Y^Mi P (p)M}, (14) 

i 

Mi = [ dad/3Ci(a,l3,t)e m e- iaQ . (15) 



In this case, if the error superoperator [acting on the 
qubit density matrix as in Eq. (|14|) ] has support in oper- 
ators that can be expanded in terms of small shifts the 
error correction will be effective. 

This is indeed the case for some classes of physical 
errors. Amplitude damping of harmonic oscillator and a 
class of unitary errors like over- rotation or under-rotation 
have been described in Ref. 0. Such errors can occur 
easily from delay or advance in the pulse application, 
for instance during the single-qubit manipulations, where 
we need to keep track of the time evolution. However, 
the physical sources of errors can be different in quan- 
tum optical and superconducting systems. While in op- 
tics imperfections in operations lead mainly to amplitude 
damping (for instance absorption of light on the beam 
splitters and phase shifters), superconducting elements 
are dephased by the Josephson junctions. So, to make 
our discussion more explicit, we consider (apart from the 
models presented already in Ref. IB) also the phase damp- 
ing of the LC oscillator. 

Specifically, let us consider the full Hamiltonian of the 
rf-SQUID, where the external flux (& x and $ x ) may be 
slightly fluctuating 



H = -2E V j cos 



7T$ X + 7T<5$ X (*) 



COS (7T$) 



TTLOq 



(0>-D x -^(t)) 2 Q 2 
2 2 



where (5$ x (t) and 8<& x {t) are the fluctuations. If we con- 
sider the oscillatory regime, i.e., $ x = ±1, we can expand 
the Hamiltonian around any integer value of <f> = k, 



H 



i 



-T7LUQ 



2 v ' V 2 
($-$ z ) 2 Q 2 
2 2 



+($-$ x )5$ x (*) 



This Hamiltonian contains terms that describe coupling 
of the flux $ to external flux fluctuations. Coupling of 
linear terms in $ corresponds to the amplitude damping 
of the oscillator, while coupling of quadratic terms in <f> 
is associated with the phase damping. (Errors generated 
by quadratic terms in the variables can be also induced 
during the single-qubit operations, as discussed in Sec. lIIII 
- they would also dephase the system.) 

Let us start with the amplitude damping of the LC 
oscillator described by the following master equation: 



t 1 t 1 t 

apa' — -a'ap— —pa' a 
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Here, 



1 



« T = ^=(*-*Q). 

For short time intervals dt, we may write 
p{t + dt) = (VTdia) p(t) [Vfdta 1 



I - ^a f adt\ p{t) (i - ^a^adt ) .(16) 



If we now compare Eqs. (fT5)l and (fTf?)) . we find the Kraus 
operators to be 

p 

Mi = VTdta, M 2 = I- -a)adt. 

The amplitude corresponding to the first operator is^ 

d(a,/3,dt) = 
i 
2 

i r 

+ 2 



6(a)<5(/3 - y/Tdt/2) ~ 6(a)6(P + y/Tdt/2) 



S(a - y/Ydb/2)8{p) - 5(a + y/Tdt/2)8[fl) 
The second amplitude C 2 is found from 



dad(iC 2 {a, (3, dt) 



J/3<S> -iaQ _ 



Tdt 



a 1 a, 



by applying inverse Fourier transform to both sides of 
this equation, which leads to 



C 2 (a,j3,dt) 



1 H 
Tdt 



™) S(a)S(P) 

^S(a)S(P)+6(a)^S(0) 



Clearly, after short time intervals, the state can be only 
slightly shifted from the code subspace. 

The amplitude damping of harmonic oscillator at fi- 
nite temperatures leads to the thermal-state solution, i.e., 
p(oo) = e~@ H . The thermal-state density matrix has 
all off-diagonal elements equal to 0. So, the amplitude 
damping gives rise indirectly also to dephasing. How- 
ever, as already noted, in Josephson-junction systems it 
is worth analyzing the pure dephasing mechanism inde- 
pendently. 

The latter process for the oscillator is described by the 
master equation 
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a'apa^a ~ —(c^a) 2 p — -p(a^a) 2 



The same procedure like for amplitude damping gives 

Mi = VTdia 1 a w — ( e 
2i V 



giy/Tdta^ a _ ^—is/Vdta* a 



Mo 



I -^{cJafdt. 



The first operator, Mi, is a sum of under-rotation and 
over-rotation of the oscillator and has been discussed in 
Ref. 0. Mi is characterized by the amplitude 



C 2 (a,(3,dt) 
Tdt 



8 J 



S(a)S((3) 



8 d 



d 2 



"' 4 S(a)5({3)+6(a)^-5({3) 



da 2 



S(a)6(f3) 



d 2 

da 2 



which for small dt should be still sufficiently localized 
around the point (0,0) to enable successful error correc- 
tion. 

Regardless of the analysis made, we could also argue 
that the amplitude C in Eq. (| 13[) is, in some sense, a 
wave function. Initially, when the state is in the code 
subspace, C(a,(3,t = 0) = S(a)5((3). Any kind of qubit- 
environment interaction generates evolution of the func- 
tion and, if the process is physical, the wave function will 
smoothly spread in time. Thus, whatever the sources of 
errors, if the state's destruction is not too fast, we should 
be able reverse the effects of decoherence. 

In the considered system, we may also reverse some 
effects of imperfections in manipulating the system. For 
instance, imprecise time measurement is equivalent to 
over-rotation or under-rotation, and imprecisely adjusted 
external flux & x causes diffusion of the state in the 
plane. 

Error recovery for this class of codes requires two steps 
corresponding to two independent errors: shifts in <I> and 
Q. The error syndrome (length of the shift) is measured 
on an ancillary qubit. First, we prepare the ancilla in the 
state 



|o) + |T)= £ |$ = .s} 



S — — QO 



If the state of the first qubit is shifted by a, we have the 
initial state of the data qubit and the ancilla in the form 

|^T^(|0) + |1». 

Then, the SUM operation is performed. Since the gate 
acts on the states \j) = J2 S |3> = 2s + j) like 



SUM:|J)|fc>-|j>li©*>> 
after the SUM operation, we arrive at 



\x+~a)(\x + a) + \x + l + a)) 




s + x + a) 



The state of the ancilla is invariant under translation by 
1, and we can omit x as it can be only or 1. The state 
of the ancilla contains the information about the shift 
only and not about the actual state of the qubit. By 
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measuring (destructively) the state of the ancilla system, 
we read out a mod 1. If we then shift the state of the 
qubit by a, the error is corrected. The same procedure on 
the Fourier-transformed state (physically the operation 
needs to be delayed by tt/2ojq - after the time roles of $ 
and Q are interchanged) yields the shift in Q that is to 
be corrected. 



B. Requirements to the circuit parameters 

The ingredients required to perform quantum compu- 
tation and quantum-error-correction routines have been 
described here for a system which may be manipulated 
with great accuracy. However, in real physical systems, it 
is rarely the case. Before we conclude our discussion, let 
us itemize the main potential problems in manipulating 
the system. 

First of all, we should briefly comment on the code- 
words. The states [Eq. ©] are clearly not physical as 
they contain components with infinite energy and are 
non-normalizable. We need to replace them with ap- 
proximate, physical codewords. In Ref. HI, the physical 
codewords are defined as superpositions of narrow Gaus- 
sians weighted by a Gaussian envelope. If the width of 
the narrow peaks [which substitute the delta functions in 
Eq. ([3])] is k and we want the states to be equally local- 
ized in $ and Q, the width of the envelope should be 
In other words, the basis of the (realistic) code subspace 
can be defined as 

OO p 

\j)oc e-^ K ^ 2s+ ^/ 2 d<S>e-^- 2s -rf/* 2 \<S>), 

s=-oo ^ 

where j € {0, 1}. There is a finite overlap of the approx- 
imate codewords |0) and |1), and the physical states can 
be confused during the logical-state readout. However, 
for a moderate expectation value of the number operator 
(a^a) w 16 (a is the annihilation operator of the oscil- 
lator), the error probability is as low as 10 -6 (the esti- 
mation of the error probability due to this overlap has 
been discussed in Ref. If we require that the mean 
energy of the states is an order of magnitude smaller than 
the energy gap A (to eliminate dephasing caused by el- 
ementary excitations), we obtain the upper limit of the 
oscillator frequency roughly 1 GHz. 

The most challenging in the realization part of the pre- 
sented procedure may be encoding. The measurement is 
assumed to be close to perfect: its duration should be 
much shorter than the characteristic time scale of the os- 
cillator r m <C 1/u)q. With present-day technology, mag- 
netometers with a sensitivity of 10 -5 $o/H z1 ^ 2 are avail- 
able. If we want to obtain approximate codewords in the 
encoding procedure, we should limit the oscillator fre- 
quency even further, down to a few MHz, which makes 
the qubit quite vulnerable to external noise. (The dimen- 
sions of the system need to be increased. Moreover, as 



the oscillation period gets longer the errors accumulated 
during each oscillation become more severe.) 

The next difficulty that appears here is the absolute 
time measurement. Once we have prepared the encoded 
state of the qubit, we need to keep track of the evolu- 
tion with accuracy determined by 1/uq. This problem 
can be overcome for a reasonably small number of op- 
erations: even if the oscillator is slightly advanced (or 
delayed), we will reset the phase of the oscillator using 
error-correcting schemes. However, if the number of op- 
erations between the error-correcting subroutines is big, 
even small imperfections may shift the phase so that ef- 
fective error recovery will be hardly achievable. 

Finally, the assumed tiny peak separation in the en- 
coded states (of the order of <l>o), resulting in the con- 
straint [Eq. J7J] might require very hardly achievable val- 
ues of the capacitance and inductance of the LC oscil- 
lator. If we multiply the separation in flux by any odd 
integer, we arrive at much weaker conditions. The log- 
ical states |0) and |1) are still peaked at maxima and 
minima, respectively, and the transformation [Eq. (|9a[) ] 
is still feasible. However, increased separation in $ re- 
sults in decreased separation of the peaks in the conjugate 
variable for which the second single-qubit gate [Eq. (|9b|) ] 
does not work as intended. To overcome this difficulty, 
it would be necessary to combine the single-qubit gates 
with squeezing. 

Certainly, if quantum computation is to be of practical 
meaning, after the calculations are over we need to read 
out the state of the register. Since the logical state is 
encoded in the flux states, using the same dc-SQUID as 
for the purpose of encoding, we may determine whether 
the state is closer to or 1 (in other words, if the flux is 
closer to even or odd value). The potential problems in 
this procedure are of the same kind like those discussed 
above: we need to perform the measurement when the 
phase of the oscillator ui^t is very close to (mod 7r) or 
7r/2 (mod 7r) for the Fourier-transformed state. Also in 
this case, the time of the measurement should be much 
shorter than 



VI. CONCLUDING COMMENTS 

We described a qubit encoded in an infinite- 
dimensional system, the rf-SQUID. Utilizing the schemes 
for quantum error correction described in Ref. [H, we 
showed how to use the entire Hilbert space instead of us- 
ing simply two states from the spectrum (as it has been 
proposed previously for the flux qubits), in this way en- 
abling error-correcting routines with a single rf-SQUID. 
In principle, one can prepare the logical state and manip- 
ulate the system in a coherent manner so that universal 
computation and error correction are possible. The set of 
gates consists of symplectic operations, for which amplifi- 
cation of errors can be avoided^ and non-symplectic that 
can be realized in the system in much more natural way 
than in quantum-optical setting. We also discussed the 
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physical errors and possible difficulties in experimental 
realizations. From the latter, we see which procedures 
should be further optimized: the encoding accuracy, if 
independent of the time necessary to perform projective 
measurement, would not lead to very strong constrains 
on the characteristic time scale of the system. 

Making use of the procedures presented here, we may 
consider also different schemes for quantum error cor- 
rection in superconducting nanocircuits, which protect 
a state of a wave-packet in an entangled state of many 
oscillators (Ref. [2J and references therein). Apart from 
the problem of quantum fault-tolerant computation we 
may use the discussed operations to various, not strictly 
computational schemes, like unconditional quantum tele- 



portation of the variables Q, i.e., not only of a two- 
dimensional subspace of the system but of its arbitrary 
wave function.— 



Acknowledgements 

We are grateful to A. Shelankov and A. Shnirman for 
inspiring discussions. This work is part of the focused 
research program "Quantum information processing" of 
the DFG. Y.M. was supported by the S. Kovalevskaya 
award of the Humboldt foundation, the BMBF, and the 
ZIP programme of the German government. 



1 M. A. Nielsen and I. L. Chuang, Quantum Computation 
and Quantum Information (Cambridge University Press, 
Cambridge, 2000). 

2 D. Gottesman, Phys. Rev. A 54 , 1862 (1996). 

3 A. R. Calderbank, E. M. Rains, P. W. Shor, and N. J. A. 
Sloane, Phys. Rev. Lett. 78, 405 (1997). 

4 A. R. Calderbank, E. M. Rains, P. W. Shor, and 
N. J. A. Sloane, IEEE Trans. Inf. Theory 44, 1369 (1998). 

5 D. Gottesman, A. Kitaev, and J. Preskill, Phys. Rev. A 
64, 012310 (2001). 

6 S. L. Braunstein and P. van Loock, Rev. Mod. Phys. 77, 
513 (2005). 

7 S. D. Bartlett, B. C. Sanders, B. T. H. Varcoe, and H. de 
Guise, Proceedings of IQC'01, edited by R. G. Clark (Rin- 
ton, Princeton, NJ, 2001) pp. 344-347. 

8 S. Pirandola, S. Mancini, D. Vitali, and P. Tombesi, Eu- 
rophys. Lett. 68, 323 (2004). 

9 S. D. Bartlett, H. de Guise, and B. C. Sanders, Phys. Rev. 
A 65, 052316 (2002). 

10 B. C. Sanders, S. D. Bartlett, and H. de Guise, Proceed- 
ings of the ICSSUR'01, edited by D. Han, Y. S. Kim, B. 
E. A. Saleh, A. V. Sergienko, and M. C. Teich, 2002 (un- 
published) . 

11 D. Gottesman and J. Preskill, Phys. Rev. A 63, 022309 
(2001). 

12 S. L. Braunstein and H. J. Kimble, Phys. Rev. Lett. 80, 
869 (1998). 

13 B. C. Travaglione and G. J. Milburn, Phys. Rev. A 66, 
052322 (2002). 

14 Y. Makhlin, G. Schoen, and A. Shnirman, Nature (Lon- 
don) 398, 305 (1999). 

15 Yu. A. Pashkin, T. Yamamoto, O. Astafiev, Y. Nakamura, 



D. V. Averin, and J. S. Tsai, Nature (London) 421, 823 
(2003). 

16 T. Yamamoto, Yu. A. Pashkin, O. Astafiev, Y. Nakamura, 
and J. S. Tsai, Nature (London) 425, 941 (2003). 

17 J. M. Martinis, S. Nam, J. Aumentado, and C. Urbina, 
Phys. Rev. Lett. 89, 117901 (2002). 

18 K. B. Cooper, M. Steffen, R. McDermott, R. W. Sim- 
monds, S. Oh, D. A. Hite, D. P. Pappas, and J. M. Marti- 
nis, Phys. Rev. Lett. 93, 180401 (2004). 

19 M. Tinkham, Introduction to Suerconductivity, 2nd ed. 
(McGraw-Hill, New York, 1996). 

20 L.B. Ioffe, M. V. Feigel'man, A. Ioselevich, D. Ivanov, M. 
Troyer, and G. Blatter, Nature (London) 415, 503 (2002). 

21 Y. Makhlin, G. Schon, and A. Shnirman, Rev. Mod. Phys. 
73, 357 (2001). 

22 A. Barenco, C.H. Bennett, R. Cleve, D.P. DiVincenzo, N. 
Margolus, P. Shor, T. Sleator, J. A. Smolin, and H. Wein- 
furter, Phys. Rev. A 52, 3457 (1995). 

23 M. Reck, A. Zeilinger, H. J. Bernstein, and P. Bertani, 
Phys. Rev. Lett. 73, 58 (1994). 

24 S. L. Braunstein, Nature (London) 394, 47 (1998). 

25 The choice of units following from Eq.([2} makes the di- 
mensionless variables fully symmetric under exchange in 
the Hamiltonian Eq.([T]). It is thus the only one for which 
the oscillatory evolution is rotation around a circle (not 
ellipse) in the phase space. Conservation of the distance 
is crucial for the one-qubit operations and hence the con- 
straint. 

26 In quantum optics the oscillator is simply a laser beam. 
The oscillatory variables q and p correspond to different 
quadratures of the electromagnetic field 



